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Abstract 

A recent global fit to the CMSSM incorporating current constraints on supersymmetry, 
including missing transverse energy searches at the LHC, BR(i? s — > and the direct 

XENON 100 search for dark matter, favours points towards the end of the stau-neutralino 
(ti-x) coannihilation strip with relatively large m\ji and 10 ^ tan/3 ^ 40 and points in 
the H/A rapid-annihilation funnel with tan/3 ~ 50. The coannihilation points typically 
have — m x ^ 5 GeV, and a significant fraction, including the most-favoured point, has 
m fl — m x < m T . In such a case, the fi lifetime would be so long that the fi would be 
detectable as a long-lived massive charged particle that may decay inside or outside the 
apparatus. We show that CMSSM scenarios close to the tip of the coannihilation strip 
for tan (3 ^ 40 are already excluded by LHC searches for massive charged particles, and 
discuss the prospects for their detection in the CMS and ATLAS detectors via time-of-flight 
measurements, anomalous heavy ionization or decays into one or more soft charged particles. 
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1 Introduction 



In the constrained Minimal Supersymmetric extension of the Standard Model (CMSSM) (TJ 
[2j, in which the soft supersymmetry-breaking parameters are assumed to be universal at 
some input GUT scale, and .R-parity is assumed to be conserved so that the lightest su- 
persymmetric particle (LSP, presumed here to be the lightest neutralino, x) is stable [3|[4], 
there are several distinct regions in parameter space where the relic x density falls within the 
range allowed by the WMAP experiment and other astrophysical and cosmological observa- 
tions [5|. These include the x — stau coannihilation strip p], the funnel region where there 
are rapid annihilations through s-channel H/A poles [TJ[T] , and the focus-point strip [8]. 

Recent LHC searches for missing-energy events [9-12 , combined with the measurements 
of Mh 13 14 and B s — > fi + decay 15 16 as well as the direct XENON100 search for 



dark matter [17], constrain severely the CMSSM parameter space. The LHC has excluded 
much of the coannihilation strip, whilst XENON100 has put another nail in the coffin of the 
focus-point strip [18). A recent global analysis finds almost equally good fits in the rapid- 
annihilation funnel and towards the end of the coannihilation strip, both being compatible 
with the LHC measurement of and other constraints [18]. The rapid-annihilation fun- 
nel is favoured by LHC missing-energy ($t) searches and Mjy, in particular, whereas the 
coannihilation region is favoured by — 2 and B s — > 

It is well known that the f\—x mass difference decreases monotonically along the coan- 
nihilation strip towards its tip at large m x . Some distance before the tip of the strip, the 
mass difference falls below m T , so that the two-body decay f\ — > x + r becomes kinematically 



forbidden, three- and four-body decays dominate, and the f\ lifetime exceeds 10 s 19,20 
In this case, any fx produced inside an LHC detector would appear as a long-lived charged 
particle that may or may not decay before escaping. The sensitivities of the conventional 
LHC searches for missing-energy events to this portion of the CMSSM parameter space 
should be examined in both cases. The LHC experiments ATLAS and CMS are capable of 
detecting long-lived charged particles, and have published limits on their possible production 
in various scenarios in which they decay outside the apparatus, including some supersym- 
metric frameworks such as gauge-mediated supersymmetry breaking (GMSB), but not the 



CMSSM 21,22 



In this paper we first show that the f\ is long-lived in a favoured fraction of the remaining 
CMSSM parameter space, and then argue that searches for long-lived charged particles 
should be important features of a comprehensive strategy for supersymmetry searches and 
future global fits. We then compare the available limits on long-lived charged particles 
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with the other LHC searches included in previous global fits, and argue that the end of the 
CMSSM coannihilation strip is nigh. The small portion that has not already been excluded 
could be explored by a combination of LHC searches for Ip? events and massive charged 
particles with the 8-TeV data taken during 2012. 



2 Summary of Results from a Global CMSSM Fit 



The starting point of our study is a recent global analysis 18 of the CMSSM including the 



ATLAS search for jets +$t events with 5/fb of data at 7 TeV (9), the (presumed) Higgs mass 



measurement Mh ~ 125 GeV 13 14 , a combination of the available constraints on BR(5 S — >■ 
fi + fi~) [l5] Q the latest XENON100 constraint on direct dark matter scattering [l7], and low- 
energy constraints including — 1 and precision electroweak data. As already mentioned, 
the XENON100 results disfavour the focus-point region of the CMSSM parameter space, 
the ATLAS $t results disfavour points on the coannihilation strip with small my2, and 
(relatively) large values of mj/2 are also favoured by the Higgs mass. 



The global analysis in 18 found two 'islands' of CMSSM parameter space with compara- 
ble values of the global likelihood: a low-mass region with m\n ~ 910 GeV, mo ~ 300 GeV, 
and tan/3 ~ 16, and a high-mass region with mi/ 2 ~ 1890 GeV, m ~ 1070 GeV and 
tan/3 ~ 45. The former lies in the coannihilation region and the latter is in the rapid- 
annihilation funnel region. As shown in Table 2 of [18], the low-mass CMSSM region is 
favoured by — 2, A^(b) and BR(5 S — > fi + fi~), whereas the high-mass region is favoured 
by M w , A e {SLD) and the ATLAS 5/fb# T analysis. 

The left panel of Fig. [I] shows how the different 'islands' found in 18 populate different 
ranges of (horizontal axis) and — m x (vertical axis). Here and in the right panel the 
red (blue) lines are the Ay 2 = 2.30(5.99) contours corresponding approximately to 68 (95)% 







CL contours, and the green stars mark the best-fit point found in 18 . To guide the eye, a 
green band shows where the mass difference between the fx and y is smaller than m T . We 
note that this mass difference is 0.42 GeV for the best-fit point and that, more generally, 
points in the lower-mass region with ~ 400 GeV, corresponding to the coannihilation 
strip, tend to have small mass differences, generally ^ 5 GeV at the 68% CL and often < m T . 
On the other hand, points in the higher-mass region, corresponding to the rapid-annihilation 
funnel, have mass differences that may extend above 25 GeV at the 68% CL. It is the small 



1 The results in 18 are very little affected by incorporation of the recent measurement of BK(B S 



/.i + fi ) by the LHCb Collaboration 16 , whose central value is quite close to the combination of pre- 



vious data, and very consistent with the predictions of the CMSSM and NUHM1 fits in [18 . See 



http://mastercode.web.cern.ch/mastercode/ for current Mastercode results. 
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Figure 1: Portions of the (m ?1 ,m ?1 — ui x ) plane (left panel) and the (m^ — m x ,tan/3) plane 
(right panel), displaying the best-fit point (green star) and the 68% and 95% CL contours 
(red and blue lines, respectively) found in a recent global analysis t 18l of the CMSSM. The 

are shaded green. 



bands with 



m x < m 7 



mass difference that makes the coannihilation mechanism sufficiently efficient to bring the 
relic x density down into the WMAP range, whereas larger mass differences are possible in 
the region where rapid XX annihilation occurs via the direct s-channel H/A poles. We note 
in the left panel of Fig. [I] a tendency for smaller values of — m x to be favoured as 
increases from ~ 300 GeV to ~ 600 GeV, as expected since m fl — m x — > as one approaches 
the tip of the coannihilation strip. The right panel of Fig. [T] shows how the different 'islands' 
found in [18] populate different regions of the (m^ —m x , tan (3) plane. We see that the points 
with tan /3 < 43 have — m x < 5 GeV, whereas those with tan f3 > 43 may have values of 
m fl — m x extending above 25 GeV at the 68% CL. The band with m fl — m x < m T is again 
shaded green. 

In the left panel of Fig. [2] we display the global A% 2 function as a function of — m x in 
the full CMSSM sample analyzed in 18 . We see that small values of m fl — m x are slightly 
favoured, though values up to ~ 20 GeV are allowed at the A% 2 ~ 1 level. However, if 
we restrict to the tan /3 < 43 region, as shown in the right panel of Fig. |2j we see that the 
preference for small — m x is much more marked. Specifically, we see again that the 
best-fit point has m fl — m x < m T , inside the range shaded green. 

The CMSSM results from 18 shown above include the constraints provided by results 
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Figure 2: T/ie difference in x 2 compared with the best-fit point as a function ofm^ —m x for 
the full CMSSM sample (left panel) and the coannihilation subsample with tan/3 < 43 (right 
panel), as found in a recent global analysis J8l The ranges — m x < m T are shaded 
green. 

from the ATLAS search for jets + J^t events with ~ 5/fb of data at 7 TeV [9] Q but not 
the constraints from searches for long-lived charged particles. However, the latter are the 
relevant constraints when — m x < m T and will become even more relevant for future 
CMSSM searches at the LHC, as we discuss below in more detail. 

We have made a similar analysis of the region of the NUHM1 parameter space at large 
mi/2 and small mo where chargino-neutralino-stau coannihilations bring the relic x density 



into the WMAP range, as commented in 18 . In most cases in this region, the NLSP is 
the fi, but in some cases it is the next heavier neutralino, x®, or another sparticle. The 
left (right) panel of Fig. [3] displays the x 2 function as a function of m NLS p — m x at values 
^ 20 GeV for points with a f\ NLSP in this portion of the NUHM1 parameter space. We 
see that in the f\ case (left panel of Fig. [3]) the x 2 function is relatively flat. This is because 
the appearance of coannihilations with the charginos and heavier neutralinos and/or other 
sparticles enables mf 



to be somewhat larger than along the CMSSM coannihilation 
strip. Note that at the best fit point of the full NUHM1, we have m? 1 — m x = 2.98 GeV and 
while the relic density is still determined by stau coannihilations, the mass difference is large 
enough so that the stau is not long-lived in this case. On the other hand, the x 2 function is 
a distorted parabola in the X2 NLSP case (right panel of Fig. |3j): at large m x o — m x the relic 
LSP density lies above the WMAP range, and at small m x o — m x it lies below the WMAP 



2 A recent ATLAS analysis of jets + 1£t events with ~ 5/fb of data at 8 TeV 11 quotes a lower limit 



m i/2 <^ 700 GeV in the range of mo close to the coannihilation region for tan /3 = 10, corresponding to 
<; 300 GeV, which does not exclude the best-fit CMSSM point found in [18]. 
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Figure 3: The difference in \ 2 for points in the NUHM1 coannihilation region discussed in 
a recent global analysis M, compared with the best-fit point, as a function of tunlsp ~ m x . 
The left (right) panel shows NUHM1 points with a f\ (yQ) NLSP, and the ranges where 
rriNLSP — m x < m r are shaded green. 

range. We also note that the minimum value of x 2 hi this coannihilation region, which occurs 
in the f% NLSP case, is ~ 1.8 higher than at the absolute minimum in the NUHM1. Because 
of the larger mass difference and the higher x 2 , we conclude that there is less reason than in 
the CMSSM to study the low-mass-difference region in the NUHM1. Accordingly, we do not 
discuss it further in this paper, though some of our later considerations within the CMSSM 
context may also apply to the NUHM1 and other models. 

3 Towards the Tip of the Coannihilation Strip 

As remarked above, the mass difference — m x required to bring the relic x density 
within the WMAP range decreases monotonically as m^j increases along a coannihilation 
strip, vanishing at its tip. Examples of this are displayed in Fig. |4| for tan /3 = 10 in the 
left panel and for tan/3 = 40 in the right panel, in both cases for the choices Aq = and 
2.5 mo (solid blue and red lines, respectively). Along these strips, ttlq has been chosen so 
the neutralino relic density lies within the WMAP range, using calculations with the SSARD 
code (23]. The locations of the tips of the strips vary with both tan/3 and Aq, but the 
qualitative behaviour of — m x as a function of the distance from the tip is similar for all 
of the coannihilation strips. 

Fig. [4] also displays the lower limit > 700 GeV quoted by the ATLAS Collaboration 
for the choices tan/3 = 10 and Aq = in a recent analysis using 5/fb of data from the 
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Figure 4: The mass difference Am = m? 1 — m x as a function of m\i% along the CMSSM 
coannihilation strips calculated using SSARD '231 for tan/3 = 10 (left panel) and tan/3 = 40 
(right panel), and for Aq = (blue lines) and 2.5 mo (red lines). The bands with m^ —m x < 



m T are shaded green. The lower limit on from the 8-TeV ATLAS 5/fbJfiT search at the 
LHC [IV is represented by maroon lines, and the lower limits from searches for the direct 

are shown as dashed and solid lines, 



and total production of metastable charged particles 
respectively, inside the green bands (see Section 6 for details). 



LHC running at 8 TeV [IT] . It was shown in [18] that a previous lower limit from ATLAS 
using 5/fb of 7- TeV LHC data M was largely independent of tan/3 and Aq. Combining 
the 7- TeV ATLAS data with other constraints on the CMSSM, that analysis also found 
m i/2 ~ 700 GeV at the A% 2 = 2.3 level, corresponding to the 68% CL in the (m ,m 1 / 2 ) 



plane, considering all values of tan /3 and Aq 18 



We see that the ATLAS limit intercepts the WMAP strips for tan /3 = 10, A = (2.5 m ) 
where m fl - m x ~ 2.8(2.5) GeV, and the WMAP strips for tan/3 = 40, A = 0(2.5m ) 
where — m x ~ 6.5 (9) GeV. We also show, again in green shading, the band where 
— m x < m T . The efficiencies and sensitivities of the current LHC$r search strategies 
need careful study in the limit of small m fl — m x , where in the CMSSM also the masses of 
the other sleptons are closer to m x than in generic regions of the (mo, m^) plane. We also 
note that the sensitivity of$r searches by ATLAS and CMS to points within the green band 
may be affected by the fact that events with a massive long-lived f\ in the final state have 
a different experimental signature. One should, in particular, consider the possibilities that 
the f\ may decay inside the detector into either one or three light charged particles, as well 
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as the possibility that it escapes from the detector before decaying. 

Clearly, the ATLAS $t search has not yet reached the band where — ui x < m T along 
any of these WMAP strips, although it is very close for tan j3 = 10. Equally clearly, future 
ATLAS and CMS searches using the full data set at 8 TeV and higher energies will reach 
this band sooner or later, depending on the values of tan (3 and Aq. 

How sacrosanct are the WMAP strips? CMSSM parameter sets lying above these strips 
would give too much cold dark matter, if one assumes that R parity is conserved (as in the 
CMSSM), that the lightest neutralino \ is the LSP, and that the Universe has expanded 
adiabatically since the epoch of thermal dark matter decoupling, with only dark matter 
and Standard Model particles contributing to its expansion (as in conventional Big Bang 
cosmology). Thus, discovery of the .R-conserving CMSSM with parameters measured to be 
above the WMAP strip would imply that either conventional Big Bang cosmology should 
be modified or some other sparticle is the LSP, e.g., the gravitino. Conversely, CMSSM 
parameter sets lying below these strips would give too little LSP dark matter. However, 
in this case one could postulate the presence of some other contribution to the cold dark 
matter density, e.g., axions, without having to rethink either Big Bang cosmology or the 
nature of the LSP. Therefore, parameter sets below the WMAP strips, in particular inside 
the green bands at lower masses, cannot be excluded absolutely. Hence it is relevant to 
consider the sensitivity of current constraints within the green band also at smaller m 1( / 2 
below the coannihilation strips. 

In the following sections, we first study the sensitivity of the ATLAS search for small 
m-fj — m x > m T in the neighbourhoods of the coannihilation strips, and then the sensitivities 
of the LHC-^t searches inside the green shaded bands where — m x < m T . 

4 Sensitivity of the ATLAS I$t Search in the Coanni- 
hilation Strip Region 

The stau-neutralino coannihilation strips, in addition to featuring a small mass difference 
between the fi and the Xi as seen in Fig. |4| may also feature relatively small mass differences 
between the x and the other sleptons, £r 7 Li f^R,L and f 2 , depending on the values of tan/3 
and Aq as well as mi/2- Since the amount of$r and the transverse momenta of any final- 
state leptons in events in which the P^r,l and fi j2 are produced depend on these mass 
differences, it is natural to ask whether there is any impact on the sensitivity of LHC searches 
for supersymmetry in the neighbourhood of the coannihilation strips. 

In this connection, we note that the granularity of the published results of LHC Ifir 
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searches in an (mo) m i/2) plane is typically 0(50) GeV , so any fine structure close to 

the f\ — x LSP boundary would not have been apparent. We also note that the ATLAS 
8-TeV^T search has a somewhat reduced reach in m^ for m ^ 500 GeV 111]. For these 
reasons, we have extended validations 18 of the previous ATLAS Ifix search in simulations 
based on PYTHIA and Delphes to study in more detail possible effects on the reach 
in the coannihilation strip region where m^ — m x = 0(5) GeV. We concentrate on the 
search for events wither but no identified accompanying leptons, since this has previously 
been found to be the most sensitive search at the large values of the ratio mi/2/mo near the 
coannihilation strips. 

We have checked the evolution of the sparticle spectrum for CMSSM scenarios with 



mi/2 — 700 GeV (similar to the 95% CL lower limit given by ATLAS for small mo 11 ), 
tan /3 = 10, 40 and Aq = 0, 2.5 mo, varying mo so that Am = m^ — m x increases from being 
adjacent to the green band where the f\ is long-lived to ~ 10 GeV. We do not see any notable 
variations in the sparticle spectrum, specifically none in the hierarchy of sparticle masses. 
On the other hand, the hierarchies of sparticle masses for tan /3 = 40 are rather different 
from those for tan /3 = 10, particularly in the Aq = 2.5 mo case. 



Fig. [5jdisplays the exclusion confidence levels we estimate using PYTHIA 24 and Delphes 25 
to simulate the 8-TeV ATLAS 5/fb Ifix search for a selection of scenarios with ra\M = 
700 GeV, the same values tan /3 = 10, 40 and Aq = 0, 2.5 mo and varying Am ^ 10 GeV. In 
these cases, the WMAP-compatible coannihilation strip has m^ — m x ~ 2.5 GeV to 9 GeV, 
as seen in Fig. |4} We see no significant differences, outside the statistical uncertainties in 
our analysis of the exclusion confidence levels, between simulations for different values of 
tan/3 and Aq. This is in agreement with the previous analysis in 18 but now with finer 
granularity in the relevant region of small m fl — m x |^j 

On the basis of this study, we assume for our purposes here that the 95% CL limit in the 



(mo, mi/2) plane established by ATLAS 11 on the basis of an$r analysis with 5/fb of 8-TeV 
data may be extended through the WMAP coannihilation strip range of m? 1 — m x down to 
the band close to the LSP boundary where the f\ becomes long-lived, simply by assuming 
that the results are independent of mx/i in this region. The vertical maroon lines in Fig. [4] 
demonstrate the corresponding impact of the IpT constraint on the CMSSM coannihilation 
strips discussed earlier. 

We see that a significant fraction of the remaining portion of the coannihilation strip not 
excluded at the 95% CL by the above analysis has Am < m T , in which case the f\ is long- 

3 Since this analysis is similar to the PYTHIA and Delphes validation described in [l8] , we do not discuss 
details here. 
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Figure 5: Estimates using PYTHIA and Delphes simulations of the exclusion sensitivity of 



the ATLAS 5/fbJ^T search along the line in the (m ,mi/2) found in 111 to yield 95% CL 
exclusion for tan/3 = 10 and Aq = 0. Our simulations are for tan/3 = 10 and 40 and the 



choices Aq = and 2.5 mo, and are plotted as functions of Am = m^ 
for tan /3 = 10 the W MAP -compatible coannihilation strip has Am ^ 
tan/3 = 40 it has Am ~ 6.5 to 9 GeV. 



- m x . We recall that 
3 GeV, whereas for 



lived. This conclusion is particularly strong for tan/3 = 10. The experimental sensitivity to 
mi/2 requires a separate, dedicated analysis in the low- Am regions, as we now discuss. 

When Am > m T , supersymmetric cascades are likely to include f\ — > r\ decays followed 
by r decays that are mainly hadronic. On the other hand, if Am < m T , the f\ decays would 
not occur promptly, and the f\ would not be registered as a conventional physics object. 
This change in the pattern of cascade decays is reflected in the populations of different 
experimental signatures, as shown in Fig. [6j The pie charts show the fractions of a pair of 
100,000 supersymmetric events generated in PYTHIA simulations that yield Ifix events with 
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no leptons (red segments) and with single leptons (green segments), as well as same- and 
opposite-sign dilepton events (SS, blue segments; OS, yellow segments), and events with no 
distinctive experimental signature (black segments). Comparing the simulation of the best- 
fit point (left) that has Am < m r with that of a point that has m fl increased so as to lie 
just outside the green band of small Am (right), we see that the best-fit point has a larger 
fraction of same-sign dilepton events (28.3% vs 17.5%) and a smaller fraction of zero-lepton 
$t events (20.0% vs 35.3%) Q These differences are due essentially to the disappearance of 
hadronic r decays from the supersymmetric cascades mentioned above. A priori, this raises 
the prospect that the efficiency and hence the sensitivity of the ATLAS $t searches discussed 
above would be reduced inside the green band where Am < m r . 



28.3%, SS 



37.0%, 1-lepton 




1.2%, OS 



20.0%, 0-lepton 
(5.5%, ATLAS) 



34.2%, 1-lepton 



13.5%, no search 



35.3%, 0-lepton 
(3.6%, ATLAS) 




2.0%, OS 
11.0%, no search 



17.5%, SS 



Figure 6: Pie charts displaying the fractions of 100,000 events generated using PYTHIA 
with parameters corresponding to the CMSSM best-fit point with Am < m r (left) and an 
analogous point with m ?1 increased so that Am > m r . Indicated in parentheses are the 
fractions of these events that survive the ATLAS triggers and selections for zero-lepton +IpT 
events according to simulations using Delphes. 

However, we then passed the event samples simulated using PYTHIA through Delphes to 
estimate which fractions of the generated events would pass the ATLAS trigger and event 
selection for the zero-lepton + lfi T category, assuming that the long-lived f\ would escape the 
detector as a massive charged particle before decaying. The results are shown in parentheses 
in Fig. [6] 5.5% of the simulated best-fit events survived the simulated ATLAS analysis chain, 
compared with 3.6% of the events with larger Am. This result indicates that the sensitivity 
of the ATLAS zero-lepton + Ifi? search for a low-Am point inside the green band is at least 

4 For completeness, we note that the fractions of one-lepton events are similar in the two cases (37.0% vs 
34.2%). 
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as high as for a point just outside the green band. This indicates that the (maroon) LHC 
limit in Fig. [4] may (conservatively) be extrapolated into the green band at low Am. 

However, a more careful consideration of prospective experimental signatures, and hence 
sensitivity, in the region where Am < m T requires a discussion of the fx decay lifetime and 
branching ratios, to which we now turn our attention. 

5 Lifetime and Branching Ratios for f\ Decay 

As was discussed in [19|, if Am = —m x > m T the dominant fx decay is two-body, namely 
fx — > tx, which occurs promptly with such a short lifetime that no fx track is detectable, 
as assumed above. However, if Am < m T the dominant decays are three- and four-body, so 
the fi lifetime is much longer, and it may decay either inside or outside the detector. 

We have recalculated the fx lifetime for the same supersymmetric model parameters as 
assumed in [l9j, namely = 300 GeV and a tl — Tr mixing ang le 6 T = tt/3 We display 
our result in Fig. [7] as a function of Am = m fl — m x . As one would expect, the f\ decays 
promptly with a lifetime ^ 10~ 20 s if Am > m T . On the other hand, as seen in more detail 
in the right panel of Fig. [7j when m T > Am > 1.2 GeV the fx lifetime is between 1 and 
400 ns, corresponding to a significant likelihood of observing the f\ decay inside an LHC 
detector, as we discuss below. We note in passing that, whereas the total f\ decay rate is 
very sensitive to Am (typically ~ Am 5 or more), it is much less sensitive to (oc l/m fl ), 
and hence the results in Fig. [7] are typical of the range of m 1 / 2 likely to be of interest to the 
LHC experiments in the near future. 

What would be the experimental signature of fx decay inside an LHC detector? In Fig. [8] 
we show results of our calculations of the dominant fx decay branching ratios. As expected, 
the dominant branching ratio for Am > m T is the two-body decay rf — > t~%. In the range 
m T > Am <: 0.8 GeV, there is competition among the three-body decays rf — > -k~v t x, 
7"f — > p~v T \ an d ^i ~^ a \ v TX (which were not considered in |19|), and the four-body decays 
rf e~v e v T x and ff — > yTVyVrX- At lower Am, the decay f-f — > n~v T x is dominant 
for Am ^ 0.16 GeV, and then rf — > e~v e v T x at Am ^ 0.16 GeV. A general conclusion, 
then, is that four potential signatures may be of interest to the LHC experiments, namely 
decays producing e~, , it~ (perhaps accompanied by one or more ir° mesons from p~ or 
di decays), and 7r~7r + 7r~ from decays. 

The right panel of Fig. [8] displays in more detail the dominant branching ratios in the 

5 Details of our calculation are given in the Appendix, where we also discuss the aspects of our calculation 
that differ from that of Jl9| . 
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Figure 7: The ?x lifetime calculated for m fl = 300 GeV and a t^ — tr mixing angle 8 T = 7r/3, 
as a function of Am = m^ — m x . The left panel covers the range 10 MeV < Am < 10 GeV 
where the lifetime is between ~ 10 12 and ~ 10~ 22 s, and the right panel shows in more detail 
the restricted range 1.2 GeV < Am < m T where the lifetime is between ~ 1 and ~ 400 ns. 
The vertical dashed lines correspond to the r, a\, p, tc and fi masses, indicated by the labels 
on the top of the figures. 
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Figure 8: The principal f\ branching ratios calculated for = 300 GeV and a Tl — tr 
mixing angle 6 T = 7r/3 ; as a function of Am = m fl — m x . The left panel covers the range 
100 MeV < Am < 2 GeV, and the right panel shows in more detail the restricted range 
1.2 GeV < Am < m T . The black, blue, orange, brown, yellow, and red lines are for the final 
states with t, ai(1260), p(770), it, p, and e, respectively, indicated by the labels adjacent to 
the corresponding curves. 
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mass range m T > Am > 1.2 GeV where, according to our previous discussion, the lifetime 
^ 400 ns, and the decay signatures are most likely to be of interest to the LHC experiments. 
Over most of this range of Am, the decay f-f — > p~v T X * s expected to dominate, with the 
next most important decay modes being f-f — > it~v T x, r-f — > e~u e v r x and ff —> \iTv^v T X- 
The branching ratio for the mode f-f — > a\v T \ exceeds 5% for Am ^1.5 GeV, indicating 
that the 7r~7r + 7r~ decay signature would have a branching ratio of a few % in this mass range. 

In the following, we discuss first and primarily the case where the f\ is sufficiently long- 
lived to escape from the detector, returning subsequently and more briefly to the case of f\ 
decay inside the detector. 



6 LHC Limits on Long-Lived Massive Charged Parti- 
cles 

So far, the most stringent LHC limits on the production of long-lived massive charged parti- 
cles have been published by the ATLAS Collaboration, based on 4.7/fb of data at 7 TeV in 



the centre of mass 22 . Upper limits are given on both the total production of metastable 
charged particles and on their direct production. In the former case comparison is made with 
the total cross section expected in a gauge-mediated supersymmetry-breaking (GMSB) sce- 
nario, and in the latter case with the cross section for producing directly three light slepton 
flavours, calculated in both cases using PROSPINO |26|. 

We use here f\ production cross section calculations obtained from PYTHIA. We see in 
the left panel of Fig. [9] that, as expected, the direct production cross section is essentially 
independent of both tan (3 and A within the range explored (tan = 10, 40, A = 0, 2.5 mo). 
We find that PYTHIA agrees to within ~ 20% with the PROSPINO calculation for three light 
slepton flavours used by ATLAS (shown as black crosses), with a difference due to the 
additional QCD corrections incorporated in PROSPINO. However, in the type of CMSSM 
scenario we discuss, most of the e,R, Jj,R decays do not lead to a fi, and the same is true 
for the decays of the heavier cl^l and T2. Therefore, in our case effectively only one light 
slepton flavour contributes to the direct production cross section, namely the f\ itself, rather 
than the three flavours considered by ATLAS. Thus, in our CMSSM scenario the effective 
cross section for the direct production of metastable charged particles shown in the left panel 
of Fig. [9] is a factor ~ 3 smaller than assumed by ATLAS. For this reason, the lower limit on 
from the absence of direct production is reduced from the value of 278 GeV quoted by 
ATLAS to ~ 170 GeV in the CMSSM, corresponding to m 1/2 > 400 GeV. This lower limit 
is displayed as dashed lines in the green bands in Fig. |4} 



13 



10 



10 



10 



1 




. tanp=10,A=0 

tanp=10, A =2.5 
A tanp=40, A =0 
x tanp=40, A°=2.5 
+ ATLAS Result 
ATLAS Limit 










\ 

•x 


X 






X 






X 

X 

I I I I I I I I I I I I I I I 



100 150 200 250 300 350 400 450 500 

' s , a JGeV] 



10 J E 



10* E 



10r 



10 



tanp=10, A Q =0 
tanp=10, A o =2.5 
tanp=40, A =0 
tanp=40, A =2.5 
ATLAS Limit 



• x 



x 

i 1 1 i 1 1 i 1 1 1 1 i 1 1 i i i 1 1 i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

100 150 200 250 300 350 400 450 500 

"WGeV] 



Figure 9: The direct (left panel) and total (right panel) f\ production cross sections as 
functions of m fl along the boundaries between the \ and f\ LSP regions in the CMSSM for 
tan/3 = 10, Aq = (blue circles) and 2.5m (green squares), and tan/3 = 40, A Q = (red 
triangles) and 2.5 mo (mauve crosses). Also shown as black dots are the upper limits on the 
corresponding cross sections from the ATLAS Collaboration and (in the left panel only) 
the result of the ATLAS calculation of the cross section for three light slepton flavours. 



On the other hand, as seen in the right panel of Fig. |9j the PYTHIA calculation of the 
total cross section for f\ production in the CMSSM yields values that depend significantly 
on tan/3 (and slightly on Aq). The larger production rate for tan/3 = 40 is because the f\ is 
significantly lighter than the en, fin at large tan/3, leading to increases in the cascade decay 
branching ratios producing f\ relative to those producing e#, The difference in the total 
cross section implies a difference in the lower limit on m fl that can be obtained from the 
ATLAS search: m fl ^ 270 GeV for tan/3 = 10, corresponding to m 1( / 2 ^ 630 GeV, and 
rrif x ^ 330 GeV for tan/3 = 40, corresponding to m\/2 ^ 730 GeV. 

These lower limits on m.1/2 from the upper limits on the direct and total cross sections, 
which are applicable in our CMSSM scenario for m fl — m x < m T , are indicated as dashed 
and solid lines, respectively, within the green bands in Fig. |4} We note that they complement 
the lower limits from^r searches. In particular, though the direct limit is weaker than the 
IpT limit on m^, the total production limit has comparable reach. We emphasize that the 
metastable particle searches exclude points in the green band that yield an LSP density below 
the estimate of the cold dark matter density based on the data from WMAP et al, which 
would have been allowed in the presence of some other contribution to the cold dark matter 
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density. In the cases studied, the excluded region of the green band does not quite extend 
to the strip where the neutralino LSP provides all the cold dark matter density indicated by 
WMAP et al. 

The vertical lines in the green bands in Fig. [4] extend from Am = only to 1.2 GeV, 
since our calculation shown in the left panel of Fig. [7] suggests that the f\ has a significant 
probability of decaying inside the ATLAS or CMS detector if 1.2 GeV < Am < m T , as we 
discuss below in more detail. 



7 Prospects for Searches for Long-Lived Staus 

We now discuss the possible experimental signatures of long-lived staus in the green band 
of interest described above, with a view to the perspectives for future searches. To this 
end, we focus on two sample scenarios: one with tan (3 = 10, A = and the other with 
tan/3 = 40, Aq = 2.5 mo, both with ~ 270 GeV and Am = m^ — m x < m T |^| 



Fig. 10 displays scatter plots of stau production events generated using PYTHIA in the 
(J3,T)) plane, where (3 is the f\ production velocity divided by c, and rj is its pseudo-rapidity 
at production. These plots are for the tan /3 = 10, Aq = case (left panel) and the tan /3 = 
40, Aq = 2.5 mo case (right panel), and each contains ~ 100,000 events. We recall that a 
typical experimental range is \r)\ < 2.5, which includes the great majority of the f% events 
generated. 




Figure 10: Scatter plots of stau production events in the 77) plane, where f3 is the f\ 
production velocity divided by c, and 77 is its pseudo-rapidity at production. The left and 
right panels are for tan/3 = 10, Aq = and tan/3 = 40, A = 2.5 m , respectively. In both 
cases, we use PYTHIA to generate event samples with m^ ~ 270 GeV and m^ — m x < m T . 

6 The production kinematics are essentially independent of Am over this range. 
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Fig. 11 displays the integrated rj distributions for the f\ in these two sample scenarios, 
showing how they are peaked around r\ = with only small tails beyond \r)\ = 2.5. There is 
a very slight tendency for the rj distribution of the fi to be more sharply peaked around zero 
in the tan/3 = 40, A = 2.5 m case (right panel), with an RMS spread Arj = 1.1 compared 
with 1.2 for the tan /3 = 10, Aq = case (left panel), reflecting a greater preponderance of f\ 
production via cascade decays. 




■o 0.0040 
1 0.0035 




Figure 11: The distributions in rj of the same f\ samples as in Fig. [70| in the left and right 
panels, respectively. 

Since the /3 distributions for the f\ are very similar in the scenarios we study, when 
integrated over rj, we just display the example tan /3 = 10, A Q = in the left panel of Fig. [12] 
It has (/3) = 0.76 and an RMS spread A/3 = 0.17. The right panel of Fig. 12 displays on a 



logarithmic scale the cumulative /3 distribution. We recall that many of the slower-moving 
produced fiS would exhibit somewhat enhanced ionization. In particular, we see in the right 
panel of Fig. 12 that a fraction ~ 1% of the fis would have < 0.25 and hence be likely to 



pass the MoEDAL 27 threshold for detecting heavily-ionizing charged particles. 



The r] and /3 distributions are promising for time-of-flight searches for slow-moving mas- 
sive particles in the large LHC detectors, so we examine in more detail the potential time-of- 
flight signatures of the produced fiS, analyzing separately the cases of the CMS and ATLAS 



experiments treating their respective detector geometries using Delphes. Fig. 13 displays 
scatter plots of the time T that the fiS produced inside the CMS (upper panels) and ATLAS 
(lower panels) detectors would take to pass through their respective time-of-flight systems, 
as functions of 77. As previously, the left panels are for the case tan j3 = 10, Aq = and the 
right panels for tan/3 = 40, Aq = 2.5 mo- As could be expected because of the larger size 
of the ATLAS detector, the spread of T values is significantly greater than for the CMS 
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Figure 12: The distribution in (3 (left panel) of the same tan/3 = 10, A — f\ sample as 
in Fig. 10, and the cumulative distribution in (3 on a logarithmic scale (right panel). 



detector. 

Fig. 14 displays the T distributions for the same examples (constrained to the ranges 
1 77 1 < 2.5 corresponding to the detector acceptances): CMS (upper panels) and ATLAS 
(lower panels), tan/3 = 10, A = (left panels) and tan/3 = 40, A = 2.5 mo (right panels). 
We see that the mean values of T are quite similar in the two models studied: in the case of 
CMS (T) = 28 to 29 ns, whereas in the ATLAS case (T) = 44 to 45 ns, and the RMS spreads 
are also similar: in the case of CMS AT — 20 ns, whereas in the ATLAS case AT = 29 to 
31 ns. We recall that the LHC has been running with a bunch spacing of 50 ns, and that 
there is discussion of reducing this to 25 ns in the future. Thus there would appear to be 
a possibility of confusion between the particles produced during different bunch crossings, 
which may increase in the future. 



In order to study this possibility, in Fig. [15] we plot distributions in a different quantity, 
namely the difference in the time-of-flight between a fi and a muon with the same initial 



momentum and pseudo-rapidity r\. Fig. 15 displays the differences in time-of-flight for the 
ATLAS detector (left panel) and the CMS detector (right panel). The distributions are very 
similar for the different models studied. As could be expected from the larger geometrical size 
of the ATLAS detector, the time-of-flight differences in CMS are generally smaller. About 20 
(12)% of the fis produced in ATLAS (CMS) would have a time-of-flight difference exceeding 
25 ns, the minimum LHC bunch spacing currently foreseen, whereas only very small fractions 
have time-of-flight differences exceeding 50 ns, the current LHC bunch spacing. 
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time to leave detector [ns] time to leave detector [ns] 

Figure 13: Scatter plots in the (T, rj) plane, where T is the f\ time- of- flight through the CMS 
detector (upper panels) and the ATLAS detector (lower panels), using the same fi samples 
as in Fig. [7o| in the left and right panels, respectively. 

8 Stau Decays inside the LHC Detectors 

Finally, we explore the possibility that the f\ may decay inside the CMS or ATLAS detector. 
As can be seen in Fig. [l6j substantial fractions <^ 10% of the produced fjs decay inside 
the ATLAS detector (left panel) or the CMS detector (right panel) if the lifetime ^ 250 
or 400 ns, respectively. Comparing the calculated f\ lifetime shown in Fig. [7] with the 
probability distributions shown in Fig. [16J we expect a non-negligible likelihood that a f\ 
with 1.2 GeV < Am < m r might decay inside the detector, particularly in the case of the 
ATLAS detector for which (T) is larger than for CMS. As seen in the right panel of Fig. [8j 
if 1.2 GeV < Am < m T a significant fraction of f\ decays are expected to be via the mode 
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Figure 14: The distributions in the time- of- flight T through the CMS detector (upper panels) 
and the ATLAS detector (lower panels) for the same fi samples as in Fig. [7o| in the left and 
right panels, respectively. 



f x — > a l v T Xi in which case the final state may contain three charged particles: 7r h + ti as 
well as the invisible neutrals v T and %■ 



Since Am < m T , and the typical f\ velocity is not large, as seen in Fig. 12, the charged 
fx decay products would typically have energies 0(m r ) and their detection would pose chal- 
lenges for the experiments. In the cases of decays into single charged particles, there would 
be a large mismatch between the momenta measured before and after the decay. In the cases 
of decays into three charged particles, they would each be very soft, accentuating the chal- 
lenge of detection. In either case, the f\ would not be detectable as a conventional physics 
object, and events producing fis would probably be classified as missing-energy events. On 
the other hand, within this general class of events, they would have the additional signature 
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Figure 15: The distribution in the time- of- flight difference between fis and muons produced 
with the same values of rj and momenta in ATLAS (left panel) and CMS (right panel). 




Figure 16: The fractions of produced f\s that would decay inside the ATLAS (left panel) 
and CMS (right panel) detectors, as functions of the f\ lifetime. The red (blue) curves are 
for the case tan/3 = 10, A = (tan (3 = 40, A = 2.5m ), respectively. 



of an anomalous charged track that changes its nature within the detector. This property 
might enable such events to be distinguished from backgrounder events. We do not attempt 
here to simulate the cases where metastable particles decay inside the detector because both 
the ATLAS and CMS searches for metastable particles employ modified tracking and object 
reconstruction algorithms whose accurate modelling lies beyond the scope of Delphes, and 
hence this work. 



As discussed earlier, the best fit found in a global analysis of the CMSSM in 18 has Am < 
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m T and hence a f\ that decays late. Fig. [17] shows a scatter plot of 10,000 CMSSM points 
with Am < m T populating bins in the (mo,mi/ 2 ) plane in the ranges < m < 2000 GeV, 
500 GeV < mi/2 < 2500 GeV. They were chosen from among the points analyzed in 18 so 



as to minimize x 2 m each bin Q The points are colour-coded according to the corresponding 
f\ lifetimes, with shorter (longer) lifetimes displayed in darker (lighter) shades. The selected 
points have f\ lifetimes in the range ~ (1, 1000) ns and hence, according to the analysis 
above, are expected either to decay with visible tracks or to pass through the LHC detectors 
without decaying. We note, in particular, that the global best-fit point has a f\ lifetime of 
0.0026 s. We also note that larger values of mi/2 at fixed mo correspond, in general, to longer 
fx lifetimes, reflecting the fact that these points lie closer to the tip of the coannihilation strip, 
where Am is generally reduced. We also note that the f\ lifetime tends to become shorter 
as mo is increased at fixed m^. This reflects the fact that these points tend to have larger 
values of tan /3, and the coannihilation strip moves away from the f\ LSP boundary as tan /? 
increases. Fig. 17 also displays A\ 2 = 2.30 and 5.99 contours (red and blue, respectively), 
corresponding approximately to the 68 and 95% CLs. We see that the points within the 68% 
CL region include both some with lower masses and shorter lifetimes and some with higher 
masses and longer lifetimes. Thus, this figure reinforces the message that a comprehensive 
strategy for exploring the tip of the coannihilation strip will need to combine searches for 
late decays inside the detector with searches for massive charged particles that leave the 
detector without decaying. 



9 Summary 

We have shown in this paper that (i) one of the favoured regions of the CMSSM parameter 
space remaining after LHC searches for$r events is towards the tip of the \— f\ coannhiliation 
strip, where the mass difference Am = m fl — m x = 0(m T ) (indeed, the best fit found in [18] 
has Am < 1 GeV), (ii) the fx lifetime typically exceeds 1 ns if Am < m T (and exceeds 
400 ns if Am < 1.2 GeV), (iii) if the lifetime is between 1 and 400 ns the f\ may well decay 
inside the detector, whereas it is likely to escape the detector as a massive, slow-moving 
charged particle if it has a longer lifetime (which might be detectable via measurements of 
its time-of-flight or anomalous heavy ionization). In the CMSSM these properties of the 
fi are largely independent of other parameters such as tan/3 and Aq, whereas a long-lived 
massive NLSP would not be as favoured in the NUHM1. 

The ATLAS $t searches used in this analysis are based on 5/fb of luminosity at 7 and 

7 The speckle and irregularities in this figure are due to the limited sample size explored in 118] . 
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Figure 17: Scatter plot of 10,000 CMSSM points with Am < m T populating bins in the 
(mo,mi/ 2 ) plane in the ranges < m < 2000 GeV, 500 GeV < mi/2 < 2500 GeV. They 
were chosen from among the points analyzed in [181 so as to minimize x 2 m each bin. These 
points are colour- coded according to the corresponding f\ lifetimes, with shorter (longer) 
lifetimes displayed in darker (lighter) brown. Also shown are Ax 2 = 2.30(5.99) contours, 
corresponding approximately to the 68 and 95% CLs (red and blue, respectively). 



8 TeV in the centre-of-mass. Both ATLAS and CMS have now accumulated > 20/fb of 
luminosity at 8 TeV, and these samples will enable them to extend significantly the reaches 
of such^r searches, reaching towards the tip of the coannihilation strip. However, the results 
summarized in the previous paragraph imply that complete coverage of the coannhiliation 
strip will require combining Jf,^ searches with searches for massive charged particles. So far, 
the published results of such searches use only 7- TeV data, so there is considerable scope for 
improving their sensitivity, and our analysis suggests this would be an interesting priority 
for ATLAS and CMS. 

However, it also follows from points (ii) and (iii) above that the standard searches fortf? 
and massive charged particles could usefully be complemented by searches for decays into one 
or more soft charged particles inside the detector, so we advocate optimizing searches for such 
decays using the 7- and 8- TeV LHC data. We present in the Appendix detailed calculations 
of the most important f\ decay modes. Simulating searches for such events would require 
interfacing PYTHIA with a f\ decay code such as that described in the Appendix. Since 
the ATLAS and CMS searches for long-lived particles use modified tracking and object 
reconstruction algorithms not included in Delphes, we have not attempted to simulate the 
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cases where such long-lived particles decay inside the detector with a lifetime between ~ 1 
and ~ 100 ns. We limit ourselves here to highlighting the importance of such experimental 
searches for driving a silver stake through the coannihilation strip of the CMSSM. 

Beyond the specific CMSSM scenario that motivated this analysis, we think that this 
study exemplifies the importance of an integrated approach to supersymmetry searches that 
combines 'traditional' Ifix searches with searches for more 'exotic' signatures. 
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Appendix: Metastable Stau Decay Modes and Lifetime 

We first discuss the f i — x — T couplings, taking into account the facts that the f\ is a mixture 
of tl and tr, and that the lightest neutralino \ is a mixture of B, W 3 and components. 
It is convenient to discuss the effects of ft — Tr mixing by referring to the couplings of the 
chiral supermultiplet partners of the two-component left-handed spinors tl and r£. In the 
MSSM, the interaction Lagrangian which contributes to the f\ — x ~ T couplings is, using 



the convention of 28 and 29 



1 rfiW — 
C = ~2dAdA A + ^9 {a) A*T^ A<«> a + h.c. , (1) 

where W is the superpotential, A{ and ipi denote the scalar and the fermionic components of 
the chiral MSSM superfields, respectively, A^ a is the gaugino field, g {a) is the gauge coupling 
constant, and T^ a is the (Hermitian) gauge group generator. The indices i and a label the 
chiral and gauge multiplets, respectively, and the index (a) labels the gauge group. In the 
context of the CMSSM, the generation mixings are not considered for sfermion fields, and 
the left-right mixing needs to be taken into account only for the third generation of sfermion 
fields. Therefore, the relevant superpotential is 

W = h T r R r L H\ , (2) 
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where h T is the r Yukawa coupling, Hi is the upper component of one of the two Higgs 
doublets 

Hi = I J J , #2 = I Tji ) 3 (3) 



and we denote their vacuum expectation values by {H\) = ui , (iff) = t>2 , (iff) = (H\) = 0. 
From eqs. ([!]) and pi), we derive the f\ — \ ~ T Lagrangian in two-component notation, as 

A nt = -K (f* R r L + f L r c L ) i/>° Hl + -j= {g'y L tiT L >! - g f* L r L \ 3 + g 'y Tl f R T c L X) + h.c. , (4) 

where ip° Hl is the superpartner of the H\ field, g' and g are the U(l) y hypercharge coupling 
and the SU(2)l weak coupling, respectively, A' and A 3 are the gaugino fields corresponding 
to the U(l) y generators, \y^ = — | or \y T i = 1, and the third component of the SU(2) L 
generator, \o 3 , respectively We note that since y T c corresponds to an antiparticle, it has a 
positive sign. 

To transform eq. Q into a four-component expression, we first define the four-component 
Dirac spinor 

r^fcY (5) 



and the Majorana spinners 



Then by using the formula 

^iPl^2 = , (7) 
where the \I/'s are general four-component Dirac or Majorana spinors, defined as 

* = I ^ or * - ' ^ 

we get 



A„t= - 9m T (f* RHl P L T + f L TP L Hi 
yzmw cos p v 

- (- tan 6 w r L BP L r - f* L W 3 P L T + 2 tan 6 w f R fP L B 
+ h.c, (9) 

where we have used h T = m T /vi, \g 2 {v\ + f|) = m^, tan/3 = v%jv\, and g' = gtanOw- 
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Finally, we rewrite eq. (|9| in terms of mass eigenstate fields. Using 

P L H 1 = N* 3 P L x°, P L B = N* 1 P L xl P L W 3 = N* 2 P L xl 

P R H l = N l3 P R x° z , P R B = N ll P R x° l , P R W 3 = N i2 P R x i, (10) 

and the relation between the mass eigenstates, fi^, and the interaction eigenstates, tl,r, 

fUflL Uf2L \ ( h 



- . , ■ , , r M . ,. vll) 



we get the T\ — x~ r Lagrangian 



'-T\-X-T 'lAl 



V2m w cos p V2 V2 

97JlT m iL N i3 - V2gttm6 w U; iR N t l ) P R 
mw cos p J 



t + h.c. 



\/2m w 

= fZxKc L P L + c R P R )T + h. C ., (12) 

where we assume the z-th neutralino is the lightest one, namely x- 

We now discuss the ?i lifetime and branching ratios, assuming that f\ is the NLSP and the 
lightest neutralino x is the LSP. When the mass difference between f\ and Xi Am = m fi ~ m xi 
is larger than m T , the two-body decay, f\ — > x T i is kinematically allowed and it dominates 
the f\ total decay rate. The two-body decay rate is given explicitly by 

IVbody = 1Q ^ m 3 { m fi +m x + m 4 T- 2m 2 f m 2 x - 2m\m 2 T - 2m 2 x m 2 T ) 5 



x 



cl| 2 + \c R \ 2 ) ((Am) 2 + 2(Am)m x - m 2 T ) - 2 (c L c* R + c R c* L ) m T m x ] (13) 



We are interested in scenarios with a metastable fi having a lifetime in the range ~ 1 ns 
to ~ 400 ns or more. When it is kinematically accessible, the two-body decay mode gives 
a lifetime many orders of magnitude smaller than 1 ns, as one can see from, for example, 
the left panel of Fig. [7j where we take the same parameters as in Fig. 2 of [19], namely a 
pure bino-like neutralino with m x = 300 GeV, Ufn = 1/2 and Ufi R = a/3/2. Therefore we 
neglect the contributions from all other decay modes when Am > m T . 



It was emphasized in 19 that the f\ would be metastable if Am < m T , so that the 
two-body decay is not kinematically allowed. We consider here the following dominant 
three- and four-body decay modes, namely, f\ — > xz/ r ai(1260), ft — » xi/ T p(770), fi — > x u t^, 
T\ — > X u t^ii^ and fi —> x v TV e z (the 'bars' for antiparticles are suppressed). These modes 
close in sequence as Am decreases. The three-body decay mode f\ —> x v r^ and the four- 



body decay modes were considered in 19 , but not the decay modes involving p and cti 
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fx j fx f fx I 

I I I 

Figure 18: Feynman diagrams for two-, three- and four-body fx decay modes. These diagrams 



are created using JaxoDraw '30l 



mesons. Just as T(t — > u T p(770)) and T(r — > z/ T ai(1260)) are larger than T(r — > v T -n\ so 
also fx — > xz/ r ai(1260) and fx —> xz/ r p(770) may be important compared to fx —> xv T Tt if 
these modes are kinematically allowed, so we include these modes in our calculations. The 



relevant Feynman diagrams are shown in Fig. 18 We have ignored the diagrams in which 
some supersymmetric particle, e.g., chargino, is in the propagator, because the contributions 
from these diagrams are suppressed by the large propagator masses compared to the diagrams 
with only Standard Model particles in the propagators. 

For fx — > X^t^i the coupling between the 7r and the W T boson is given by chiral 
perturbation theory (see, for example, |31|), as 

Cw-, = - 9 -^ w; d^- + h.c. , (14) 

where f^ ~ 92.4 MeV is the pion decay constant and V u d is the quark CKM matrix element. 



Apart from the coupling given in eq. (12), the only other coupling involved in this mode is 



the W-lepton-neutrino coupling given by the Standard Model, 



Cw-i-v = -^=Wp^P L l~ +h.c. (15) 



The decay rate for this process is 

£((Am; 
1287r 3 m 



Tl 



1 

21 2 



(qj-miy 



(q*-miy + {m T r T y Qf 
x [(\c L \ 2 q 2 f + \c R \ 2 m 2 T ) ((Am) 2 + 2m x Am - q 2 ) - 2 (c L c* R + c R c* L ) m x m T q 2 ] , 

(16) 

where Gf is the Fermi constant. In doing the phase space integration, we have employed 
the same substitution as in [19] , namely 

q) = (Am) 2 - ((Am) 2 - mf) x, (17) 
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where the index / denotes the massive particle in the final state other than the Xi namely, 
/ = ii here, and it denotes oi, p, n and e, respectively, in the other three-body and four-body 
decay rates. 

Forfi -> x^ai(1260) and n -> x^p(770), we use the ^-0^(1260) and W ± -p T {770) 



couplings suggested by the idea of meson dominance 32 



C = (^'" - <^ar M ) + t.c. , (18) 

where g p is the pixix coupling constant, and w p and w ai are phenomenological parameters. 



There is one difference between our treatment of w p and the way done in 32 : we take the 
id n as a phenomenological parameter to be fixed by the r decay data, while in [32] w p — 1 
exactly. This is equivalent to saying that we introduce one more phenomenological parameter 
as compared to [32] , as we get the value of w p by equating the experimental decay rate data 
for r — > ix~ix Q v T given in k33j (r 14 ) to the tree-level two-body decay rate of r — > v T p{770), 



r<2 T/2 „„2._,2 4 
r _ ^F V ud W p m p m w . 2 2 \2/ 2 ,r. 2\ / 1 qN 

Similarly, we get the value of w ai by equating the sum of the decay rates of r — > -k~2ii v t 
and r — > 7r~7r + 7r~z/ r given in [33] (r 20 and r 62 ) to the tree-level two-body decay rate of 



r — > u T ai (1260), which is obtained by the substitutions m p — > m ai and w p — > w ai in eq. (19). 
The decay rate expression for f\ —> xi/ T p(770) is 

1^^x^(770) - — 7^2 — / ax [{{Am) - q f ) {{Am) + 2m x ) - q f \ 

IzSir^mi qi [mi — m„ r ) Jo 



1287r 3 m| i 5(2 ( m 2 _ m 2^y 

1 (q) -m 2 p ) 2 {q 2 f + 2m 2 p ) 

X 



(qf - m%y + (m r r r ) 2 9/ 
x [(I c l| 2 <?/ + l c i?| 2m r) ((Am) 2 + 2m x Am - gj) - 2 (c L c^ + c R c* L ) m x m T q 2 f ] . 

(20) 

By substituting m p — > m ai and w p — > w ai in the above equation, we get the decay rate for 
fi -> xz/ r ai(1260). 

In calculating the four-body decay rates f x — >■ xv T Vpii and fx — > xz/ T z/ e e, we take the 
leading-order result in the 1/m-w expansion, since we are interested in scenarios where Am <C 
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mw- The result is 



4— body 



G\ ((Am) 2 - m\ 



96(2vr) 5 m 



x 



x 



dx [((A 

12mfqj log 



qj) ((Am) + 2m x ) 2 - -j 

9 / 



2 

mf 



™? ) (?/ - 8m 2 gJ + mf ) 



(g 2 - m 2 ) + (m T r T 

c L \ 2 q 2 f + \c R \ 2 m 2 T ) ((Am) 2 + 2m x Am - qfj - 2 (c L 4 + c fl c^ 



(21) 

where m; = m e and m; = m^ for e and /i final states, respectively. 

The branching ratios for all decay modes are shown in Fig. [8j where the same parameters 
as for Fig. 2 of [19] are used. 

Apart from the inclusion of the decay modes fi — > x z/ rOi(1260) and f\ — > x z/ rp(770), we 
found a couple of differences between our results and those given in |l9] , in which a bino-like 
neutralino is considered. In terms of the notation there, we get a different expression for 
from eq. (3) of 1 19] , and we note that our three-body decay rate expression Yf^^^ given 
m eq. (|16[ is different from eq. (A3) of [l9) Q 

Our results for the f\ lifetime and decay branching ratios are compared with those of [19] 
in Fig. 19 The effect of the difference in compared with [19] is limited to the interference 



term oc m T , which is important mainly when Am is small and 0(m T ), as seen in Fig. 19 
The difference in the expression for Tf^^^ is responsible for the difference visible in Fig. 19 



when Am ~ 200 MeV. The contributions from the three-body decays T\ — > xz/ T ai(1260) and 
Ti — > X^i-p(770) are important for Am ^ 1 GeV, with a combined branching fraction ~ 50% 
when Am ~ 1.75 GeV, and still > 30% when Am ~ 1.25 GeV. This feature is also visible 
in Fig. [T9 
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